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' We examine the neutralino dark matter (DM) phenomenology in supersymmetric scenarios 

with nonuniversal Higgs masses (NUHM) at the gauge coupling unification scale that can 

r^^ • acommodate a light Higgs boson, where the correct relic density is obtained mostly through 

t*H| the annihilation into a pseudoscalar A. Our analysis shows that most part of the A pole region 

can produce detectable gamma-ray and antiproton signals. We further focus on uncertainties 

flj ' influencing the results in indirect and mainly direct detection. 



Introduction 

>\ 
CN , 1.1 The model 

a^ 

^^ ' One of the major experimental constraints on the constrained-MSSM parameter space comes 

from the LEP-2 Hmits on the Hghtest higgs boson mass. In particular, LEP-2 set a lower bound 
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of 114.4 GeV for this mass'"*^, excluding the largest part of the model's viable parameter space. 
However, strictly speaking, this bound only applies to the SM. It comes mostly from searches 
in the Higgsstrahlung channel, which in the case of the MSSM is actuallynot identical to the 



SM one. In particular, aMssuie'^e — )■ hZ) = sin^(/3 — a)o"sM(e^e — > hZy^^. Hence, the LEP2 
bound rrifi > 114 GeV applies to the MSSM only if sin^(/3 — a) = 0(1). This is the case in 



;J| ■ cMSSM/mSUGRA scenarios. 



This limit can actually be partially circu mvent ed once some of the cMSSM constraints are 
relaxed. In particular, it has been pointed ouv ^^^ that relaxing the requirement for higgs mass 
universality at the GUT scale can effectively reduce the sin^(/3 — a) factor, leading to a smaller 
cross-section and thus to weaker bounds on the lightest higgs mass. 
Our particular model^^ is characterised by the following parameters 

mi/2, Aq, sign(/i), tan/3, mo, mjj^, mjj^ (1) 

where the GUT-scale common scalar mass thq concerns all scalars but the two Higgs bosons. 

We examined the neutralino dark matter - related phenomenology of this model, notably 
the behavior of the relic density over some region of the parameter space, the prospects for 
indirect detection as well as the constraints coming from direct DM detection experiments. 
Viable parameter space points have to satisfy a number of constraints: 

- Higgs boson mass limit: In the non decoupling region where the A boson becomes very 
light the lower limit of nih goes down to 93 GeV or even lower. We consider that the parameter 
space with sin^(/3 — a) < 0.3 (or, sin(/3 — a) < 0.6), and 93 < mh < 114 is in agreement 
with the LEP2 limit . Consequently, the coupling of the heavier Higgs boson to the Z boson 



{gzzH oc cos(/3 — a)) becomes dominant and this makes the heavier Higgs boson SM - like, so 
the LEP-2 114 GeV hmit starts applying for the heavier CP-even higgs boson. On the other 
hand, in the decoupling region sin(/3 — a) ~ 1, which means that the 114 GeV limit applies to 
the lightest higgs. Given the fact that there exists an uncertainty of about 3 GeV in computing 
the mass of the light Higgs bosonl^, we acce pt a lower limit of 111 GeV. 

- Br{b -^ s-f) constraint: We demand'2'12! 2.77 x 10""^ < Br{b -^ sj) < 4.33 x 10"^. 

- Br{Bs — 7- /i^/i^) constraint: We further impose the important Br{Bs — > /i"^/U~) constraint 
coming from CDF, Bi{Bs — )• fi^fi^) < 5.8 x 10~^ (at 95%C.L.), which has recently been 
improved to < 4.3 x 10"*^ at 95% ClD^]. 

- WMAP constraint : In computing the relic density constraint, we consider the 3a limit 
of the WMAP data-lS 0.091 < Q.cDhih'^ < 0.128. Here Q-cDuh? is the dark matter rehc 
density in units of the critical density and h = 0.71 zt 0.026 is the Hubble constant in units of 
100 Km s~^ Mpc~^. We use the code micrOMEGAS^^ to compute the neutralino relic density. 



1.2 Dark matter detection 

There exist two main modes of dark matter detection, usually referred to as "indirect" and 
"direct". Indirect detection is based on the principle that if DM can annihilate (or decay), 
in the early universe in order to give the measured relic density, this process should also occur 
today throughout the galaxy (and beyond), so we could hope to detect its annihilation products: 
gamma-rays, positrons, antiprotons and neutrinos. Direct detection of DM relies on the fact 
that WIMPs may interact with (scatter on) ordinary matter. This scattering is an in principle 
measurable effect and indeed a huge effort is currently being developped worldwide to measure 
potential signals coming from DM scatterings upon large underground detectors. 

Concerning indirect detecion, in this work we compute the gamma-ray signals at intermediate 
galactic latitudes^^ in the spirit of elliminating as much as possible uncertainties coming from 
the DM "halo profile" (i.e. its distribution in the galaxy) as well as background contributions 
to the spectrum. These gamma-rays can be detected by the Fermi satelite^^. For the case of 
antiprotons, we compute the prospects for detecti on in th e AMS-02 mission L!^. We adopt a 
semi-analytical treatment of the diffusion equation^ " l ^ ^ ^ l, presenting results for the so-called 
MAX propagation model. 

Finally, we compare the neutralino-nucleon s piii-iiidep endent scattering cross-section to some 
of the tightest bounds available in the literatur d ^ ^ i As it has been pointed out, uncertain- 
ties are not absent in direct detection as well. These can be of a twofold nature: First of all, 
when giving exclusion bounds on the (m^,cj^(p„)) plane, a set of astrophysical assumptions (as 
well as some assumptions on the passage from the nuclear to the nucleonic level) have already 
been made. These assumptions are thought to have a small impact on t he results, the mangitude 
of which depends, among other factors, on the considered mass range^ ^ ^ ^ ' l ^o | Secondly, 
there are often uncertainties in the cross-section computations performed by theorists in specific 
models. In our case, what is of relevance is one of the parameters entering the passage from the 
partonic to the the WIMP-nucleon scattering cross-section, denoted by fxs- This parameter is 
actually related to the strange quark content of the nucleon. Its value can be either measured or 
estimated through lattice QCD methods. The DarkSUSY^^^ code, which was used to compute 
the cross-section, adopts a default value of /t^ = 0.14. However, recent lattice simulations ■'■^'^ 
point towards much lower values, of the order of 0.02, being even compatible with zero. The 
effect of this uncertainty has been quantified'^ and is known to range from negligible to very 
large, depending on the specific mechanism driving the scattering cross-section. In what follows 
we shall quantify the effect of this uncertainty showing that it is really crucial in assessing the 
viability of our models. 
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Figure 1: Viable points in the 7TI1/2 — ttia plane. Neutralino masses of ~ 55 — 65 GeV correspond to the Light 
Higgs Boson region. Detectability of the photon and anti-proton signals are represented by D^ and D'^^j^^ lines. 



2 Results 



2.1 Indirect detection 



We performed'-^ two scans in the model's parameter space: 
- In the first one, we fix tan/3(= 10), tuq = 600 GeV, Aq 



-1100 GeV, sign(^) > 0. Then, 



we vary the mass parameters rrih^ (0 < rrih^ < ttiq) and mh^ ( — LSttiq < rrih^ < — O-Smg) to 
obtain hght neutrahno dark matter consistent with fight Higgs masses {niH^A ^ 250 GeV) at 
the electroweak scale. ^ and niA are derived quantities. We note that the high /x parameter 
values obtained in this scenario correspond to an essentially pure bino LSP. 
- In our second scan, we fix mo at the very similar value mo = 600 GeV, while ^q = —1000 GeV 
is chosen to make 6 — )■ S7 less restrictive. We set nih^ (= 2.4?7io) and vary rrih^ (— O.Smg < ruh^ < 
0.1) with raii2 to obtain the WMAP-compatible regions for neutralinos. Our scan renders very 
small /i values (150 < /i < 300), consequently the LSP can have large Higgsino components. 
In both cases, we look for points satisfying the relic abundance requirement while passing all 
constraints previously mentioned. The viable points are then scattered on the m-i /2 — m-A plane 
as red dots. 

Our results can be seen in figures [T] and [2j Appart from the WMAP-compliant points, 
we demonstrate in the same plots regions excluded by other constraints (gray regions), higgs 
mass isocontours, as well the parameter space regions that can be probed during a 3-year data 
acquisition period for Fermi and AMS-02: all points lying inside the two parallel lines in the 
case of figHJand on the left or below the lines in figi2]can -in principle- be probed. 

In the first set of scenarios, there are mainly two mechanisms that can generate the correct 
relic density: quasi-resonnant annihilation through a. Aox H pole, extending along the direction 
of the line where 2m ~ mA,mH, or the light Higgs pole at low mi/2 and along the tjia 
direction. 

In the second set of scenarios the neutralino self-annihilation cross-section is enhanced kine- 
matically as before (i.e. we are once again sitting near the A pole), but moreover the neutralino 
acquires a non-negligible higgsino component which enhances its couplings to the higgs bosons. 

In both scenarios, we see that the detection prospects are quite good. Significant portions 
of the viable parameter space can be probed. 

In the first case, perspectives are actually good along the A-pole. On the contrary, we see 
that the light higgs pole seems to be completely invisible in both channels. This is due to the 



> 

O 150 

< 



NUHM 

tanp=10 
in„=600 GeV, A„=-1000 GeV, 

150<^<300 (GeV) 




sin(P-a) 



50 



260 280 

m„iGeV) 



300 



Figure 2: Same as Fig[Tl except that the light Higgs boson zone is shifted to larger neutrahno mass values. 
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Figure 3: {m^o^a o_j^) combinations along with the relevant exclusion limits from direct detection experiments 

for /j^ = 0.14. Points lying above the lines are excluded according to the published limits. The light blue 
and the dark blue points represent the light Higgs boson regime for the two scenarios at stake. 



fact that whereas resonnant annihilation is an efficient mechanism in the earl y univ erse, the 
cross-section for this process tends to zero as the neutralino velocity does so^ ^ I This is 
actually the case at present times, which are of relevance for indirect detection. On the other 
hand, annihilation through a pseudoscalar is not that sensitive to changes in the WIMP velocity 
l^, so the cross-section remains relatively high even at present times. 

In the second set of scenarios, the prospects are actually even better. On the one hand, 
any interference of CP-even higgs bosons is negligible. Moreover, the fact that in this case the 
neutralino has a significant higgsino component enhances its couplings to the higgs sector, an 
effect which is practically insensitive to velocity changes. So, the cross-section remains quite 
stable at present times. 



2.2 Direct detection and associated uncertainties 

As a final step, we computed the spin-independent neutralino-nucleon scattering cross-section. 
The results can be seen if figures [3] and HI for two different values of the /t^ parameter, namely 
the default DarkSUSY value of 0.14 and a reduced value of 0.02 respectively. 
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Figure 4: As in fig[3]but for fr 
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From figi3]we see that our NUHM model seems to be hopelessly excluded, especially when 
it comes to the Light Higgs Scenarios. The situation seems to be however quite different once 
we look at figlH 

A comparison of the two figures demonstrates the importance of the uncertainty in the 
value of Jts- We see that once we reduce its value from the default one towards lower values 
the cross-section is reduced by something like an order of magnitude. We further remind that 
/Ts estimates are even compatible with zero. Moreover, once a set of astrophysical or nuclear 
uncertainties is taken into account (local density, velocity distribution, nuclear form factors), 
the experimental limits may be considered to bare an additional uncertainty of a factor 3 — 4. 
We thus see that stating whether a model is excluded by direct detection data or not may be a 
tricky issue. All potential sources of uncertainty should be examined before definitively ruling 
out a model. Fortunately, a large effort is being devoted by several groups in quantifying these 
uncertainties and incorporating them in a systematic manner, both in the limits published by 
experimental collaborations and in the calculations performed by theorists. 
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